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Abstract 

A unified equation for the non-singlet spin dependent structure function 
gi'^{x,Q'^) = g'l{x,Q'^) — (7"(x, Q^) which incorporates the complete leading or- 
der Altarelli-Parisi evolution at finite x and double logarithmic ln^(l/x) effects at 
X — > is formulated. This equation is solved assuming simple phenomenological 
parametrisation of the non-perturbative part of the structure function. Reasonable 
agreement is obtained with the recent data obtained by the Spin Muon Collabora- 
tion. Predictions for the small x behaviour of gi^{x,Q'^) in the kinematical range 
which may be relevant for the possible HERA measurements are given. The con- 
tribution to the Bjorken sum coming from the region of very small values of x is 
quantified. Extrapolation of gi{x, Q'^) to the region of small x and small is also 
discussed. 



1. Introduction 



Inelastic scattering of polarised leptons on polarised nucleons is a powerful tool to study 
the internal spin structure of the nucleon. Measurements on proton, deuteron and neu- 
tron targets allow verification of sum rules like e.g. the Bjorken sum rule, |^, which is a 
fundamental relation of QCD and which refers to the first moment of the non-singlet spin 
dependent structure function, 5f"*(x, Q^). Here the variables x and are conventionally 
defined as x = Q'^/{2pq) and = —q^ where q and p denote the four momentum trans- 
fer between the leptons and the four momentum of the proton respectively. Evaluation 
of the sum rules thus requires knowledge of spin dependent structure functions over the 
entire region of x. Since the experimentally accessible x range is limited (presently x > 
0.003, 0), extrapolations to x = and x = 1 are unavoidable. The latter is not critical, 
giving only a small contribution to the moments. However, the small x behaviour of gi{x) 
is theoretically not well established and results on the moments depend drastically on the 
assumptions made for this extrapolation, cf. Theoretical studies of the problem are 

thus awaited for. This is even more important in view of a possible future polarising the 
proton beam at HERA, [Q, and thus investigating the polarised electron-proton collisions 
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at low X which may reveal a new dynamics. 



It has recently been pointed out that the small x behaviour of both singlet and non- 
singlet spin dependent structure function gi{x,Q^) is controlled by the double logarith- 
mic terms i.e. by those terms which correspond to powers of a,ln2(l/x) The 
double logarithmic terms also appear in the non-singlet structure function Q^) = 

F2 (x, Q^) — -^^(x, Q^) but the leading small x behaviour of the F!^'^ which they generate 
is overriden by the (non-perturbative) contribution of the A2 Regge pole 0. They can 
however generate the leading small x behaviour of the spin dependent structure functions 
gi where the relevant Regge poles are those which correspond to axial vector mesons P,p!0 



expected to have low (~ 0) intercept. Let us recall that the Regge pole model gives the 
following small x behaviour of the structure functions F^ '^lxjQ'^) and g\{x,Q'^) where 
g\{x,Q'^),i = s,ns denote either singlet Q^) = gi{x,Q'^) + fi'"(a;, Q^)) or non-singlet 

{gi'^{x, Q"^) = gi{x, Q"^) — (7" (x, Q"^)) combination of structure functions: 

F2"^(a;,g2) _a;i-°A2(o) (1) 



<7i(x,Q')~x-'W (2) 

In eqs (|T]) and (0) aA2(o) ~ 1/2 is the intercept of the A2 Regge pole trajectory and ag^nsiS^) 
denote the intercept of the Regge pole trajectory corresponding to the axial vector mesons 
with / = or J = 1 respectively for which we expect as^ns{^) < 0. For ^^^(O) ~ 1/2 
the resummation of the double ln^(l/x) terms does not affect substantially the small x 
behaviour of F2(x, Q^) and its leading behaviour at low x remains to be given by equa- 
tion (|l]). On the contrary, the Regge behaviour of the spin dependent structure functions 
g\{x^ Q'^) is unstable against the resummation of the ln^(l/x) terms which generate more 
singular x dependence than that implied by eq.(|D for ^^^^^(O) < 0. 

The double logarithmic ln^(l/x) effects go beyond the standard LO (and NLO) QCD 
evolution of spin dependent parton densities |ll-|l^- They can be accomodated for in the 



QCD evolution formalism based upon the renormalisation group equations after including 
the complete "small x" anomalous dimensions with the ln^(l/a;) terms resummed to all 
orders []TB[. The purpose of this paper is to develop alternative formalism based on unin- 



tegrated distributions in analogy to the unpolarised case [|T9|,|20| but for simplicity limited 
to the non-singlet structure function Q'^). The content of the paper is as follows: in 

the next section we recall the basic integral equation for the non-singlet combination of 
the unintegrated quark distributions which resums the double logarithmic ln^(l/x) terms. 
In Sec. 3 we formulate the unified equation which incorporates the effects of the small x 
resummation and the LO Altarelli-Parisi evolution. Numerical analysis of this equation is 
presented in Sec. 4 starting from the simple parametrisation of the nonperturbative part 
of the spin structure function g"^^ . Possible extrapolation of the spin dependent structure 
function to the region of low is also discussed in this section. Finally in Section 5 we 
give a summary of our main results. 



2. Double logarithmic \t?(1/x) resummation 

The double logarithmic terms in the non-singlet part of the gi{x, Q"^) (or of the F2{x, Q"^)) 



are generated by ladder diagrams with quark (antiquark) exchange, Fig.l pTf-p 
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To be precise for the spin dependent struc- 
ture function (yfi(x, Q^) one has also to in- 
clude certain class of non-ladder diagrams 0. 
In the non-singlet case however their contri- 
bution is non-leading in the large limit 
{Nc is a number of colours). It turns out 
to be numerically small for = 3 and 

will be neglected. The contribution of non- 
ladder diagrams is however non-negligible for 
the singlet spin dependent structure func- 
tion [0. In what follows we shall con- 
sider only the non-singlet structure function 
grix,Q') = ^7f(x,Q2) _ <(x,g2). For 
notational simplicity we shall drop the su- 
perscript 'ns' i.e. we will denote Qi'^ as gi 
etc. 



x'p, 



mrrnTTiTTnT 



mrrrfTTimnnnf 



rnrrririTrrrinnf 



rrmrrrinmnnnf 



Figure 1: An example of a ladder dia- 
gram generating double logarithmic terms in 
the non-singlet spin structure function gi. 



When discussing the leading (double) logarithmic ln^(l/a;) resummation effects in the 
spin dependent structure function gi{x,Q'^) it is convenient to work with the (combina- 
tion) of the unintegrated spin dependent quark distributions, /(x',fc^): 

f{x\e) = \{Ux\e)-Ux\e)) (3) 



where x' is the longitudinal momentum fraction of the parent proton carried by a quark 
and fc^ denotes the quark transverse momentum squared. The sum of double logarthmic 
ln^(l/x) terms corresponding to ladder diagrams is generated by the following integral 
equation for the (unintegrated) structure function f{x',k'^): 



where 




In this equation only the leading part of the splitting function Pqq{z) in the z = limit is 
retained. The variables fc^(fc'^) denote the transverse momenta squared of the quarks, /cq 
is the infrared cut-off and the inhomogeneous term f^^\x', k"^) will be specified later. The 
integration limit k'^/z follows from the requirement that the quark virtuality is controlled 
by the tranverse momentum squared. The structure function gi{x,Q'^) is related in the 
following way to /(x', k"^): 

9iM') = g?ix) + /^^ —fix' = x(l + ^),P) (6) 

where gf^"^ (x) is the nonperturbative part of the structure function and 

= Q\- - 1) (7) 

X 

To be precise, the kinematic limit for the integration over dk"^ should be set equal to 
but we have checked that the results are relativly insensitive to this change. 



3 



The relation x' = x{l + k'^/Q'^) follows from the on-shell condition imposed upon the 
struck quark i.e. (g + k)'^ = where k is now the four momentum of the probed quark. 
This can be seen after we decompose this four momentum into basic light-like vectors p 
and q' = q + xp: 

k = x'p- I3q' + kt (8) 

and assume that integration over d(3 which is implicit in the definition of the distribution 
/(x', k"^) is restricted to the region (3 « 1. In the limit fc^ << we do, of course, get 
x' = X. 

The origin of the nonperturbative part Qi^^ (x) can be viewed upon as originating from 
the non-perturbative region fc^ < /cq, i.e. 

9?\x) = f^fix,e) (9) 

The origin of the double logarithmic ln^(l/a;) terms in gi{x, Q"^) can be traced to the fact 
that the conventional single logarithmic terms coming from the logarithmic integration 
over the longitudinal momentum fraction z are enhanced by the logarithmic integration 
over the transverse momentum up to the z dependent limit k"^ / z in eq.(^) and up to the 
x dependent limit = Q'^{l/x — 1) in eq. Those double logarithmic ln(l/a;) terms 
will generate the leading small x behaviour of the structure function gi{x^Q'^) provided 
that the input structure functions g'i\x^Q'^) and f^^\x^k'^) are non-singular for x — 0. 
For fixed (i.e. non-running) coupling otsik"^) — >■ <5s the small x behaviour is then found to 
be given by 

g^{x,Q^) ^ x-"^ (10) 

where 

A = 2v^ (11) 
and where we have neglected in eq.(0) the slowly varying logarithmic factors. 



3. Unified equation incorporating Altarelli-Parisi evo- 
lution and the double \v?{llx) resummation. 

In order to derive the unified equation which will include complete LO Altarelli-Parisi 
evolution and the ln^(l/x) terms resummation let us at first confront eqs (^ and (^ 
with the LO Altarelli-Parisi evolution equation for the non-singlet structure function 
gi{x,Q'^) |[Tl|, |l^. This equation can be written in the following form in terms of the 
unintegrated distribution /(x, /c^) : 

a.(P) [4 f^dz{l + z^)gf\^~2zgf\x 
3 



f{x,e\ 

as{k') 



27T 



2tt 



z 1 — z 

4 f^dz{l + z^)f{^,k^)-2zf{x,k'^) 



^ + :-ln(l-x)]gf\x) + 



3 Jx z 



+ [^ + ^Hi-x)]fix,e) 



and 



g^{x,Q') = g^^\x)+ [ 

Jk 



d^ 





f{x,e) 



(12) 
(13) 
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Equivalence of the eqs (|T^,|TB|) to the Altarelli-Parisi equation can be easily seen by ob- 
serving that the left hand side of eq.(|l2D is equal to 



k 



2dgi{x, k^) 



and its right hand side is : 



2n 



4 

3 Jx 



dz{l + z^)g,{%k^)-2zg,{x, P) 
z 1 — z 



+ t: 



^ln{l -x) ) gi{x, k'^) 



As usual the small x behaviour of the solution of the Altarelli-Parisi equation depends 
upon the small x behaviour of the input structure function. If the input structure function 
gi (x) is not singular at a; ^ then the leading small x behaviour of the structure function 
gi{x, Q"^) obtained from the solution of the Altarelli-Parisi equation is given by: 



gi{x,Q^) ~ exp[2^aQ')ln{l/x)] 



where 



4 j-Q" c/Pa,(P^ 



ftp 



(14) 



(15) 



A;2 2tt 

This leading small x behaviour is generated by the approximate form of the Altarelli- 
Parisi equation in which at z — > we keep only the leading term of the splitting function 



Pqq{z). This approximate equation has the following form 



f{x,e) = f^'\x,k')+a,{k') f 

J X 



1 dz /-fc' dk 



n 



A-2 



A;'2 



X 



k 
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(16) 



with the structure function gii^x^Q"^) given as usual by eq.(|T3D. This equation has the 
same structure as eq. (|^) except for the upper limit of integration over dk''^. It is therefore 
evident from eq. ([16D that the LO Altarelli-Parisi evolution is incomplete at low x since it 
does not generate the double but only the single ln(l/x) terms. 



Equation (|^) generates correctly the leading small x behaviour but it is inaccurate 
in describing the evolution in both for very small and, of course, for large values of 
X. In order to have the correct evolution one should thus include in the formalism 
the complete splitting function Pqq{z) and not only its leading component ai z —>■ 0. In 
order to have at the same time the possibility to generate the double logarithmic ln^(l/a;) 
terms at low x one has to modify eq.(|T2|) and change appropriately the upper limit of 
integration over transverse momentum from k'^ to k'^/z (as in eq.(^) in that part of the 
QCD evolution which is controlled by the leading part of the splitting function Pqq{z) at 
2 — >• 0. The corresponding equation then reads: 



fix',k') = f '\x',k') + asik') ['^ [[ 

Jx' Z Jki 



1 dz i-k^z dk'^ , a;. 



+ 



as{k^) dk'^ 



where 



27r A;'2 



f'\x',e) 



Osik^ 
2tx 



4 f^dz{z + z^)f{^,k^)-2zf{x',k'') fl 



3 Jx' z 



1-z 



+ (| + -Ml-x'))/(x',fc2) 



.(0)/ 



(17) 



(18) 
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and the function ds{k^) is given by eq.(^. The structure function gi{x,Q^) is related 
to /(x',A;^) by eq.(|). Equations (|1^,(|18D and (||) which treat both potentially large 
logarithms ln(l/x) and In(Q^) on equal footing will be the basis of our analysis. Similar 
treatment of the unpolarised (singlet) structure function and of gluon distributions which 
combine the Altarelli-Parisi and BFKL evolution(s) has been proposed in refs |IU|,pD 



Finally let us observe that we will need the unintegrated distribution only in the "per- 
turbative" domain, k"^ > k^, while the non-perturbative contribution is parametrised in 
terms of the input distribution gi'\x), cf. eqs (H), (^ and ([T8|). 



4. Numerical results 

We solve equation (|1^ assuming the following simple parametrisation of gf'\x): 

g?\x) = ^gA{l-x)' (19) 

where g^ is the axial vector coupling constant (we set gA = 1.257). In the small x limit we 
have gi^\x) const that corresponds to the assumption that aAi(O) = where a;^j(0) 
is the intercept of the Ai Regge pole trajectory (cf. eq. (0)). The function gi'\x) given 
by eq. ([I9|) does also satisfy the Bjorken sum rule, i.e. 

/'rfx,r(x) = f (20) 
Jo 

The parameter k^ is set k^ =1 GeV^. 

In Fig. 2 we compare the calculated gi{x,Q'^) with recent measurements of the SMC 
which extend to presently lowest values of x. In view of the simplicity of the parametri- 
sation (|T^) with no free parameters, agreement with the data can be regarded as quite 
satisfactory. Parametrisation of gi'\x) can certainly be made more flexible, with phe- 
nomenological parameters to be fitted to the data but large experimental errors make 
such refinements irrelevant at present. 

In Fig. 3 we show the structure function gi{x,Q'^) as the function of x for =10 
GeV^ obtained from eq.(^ with /(x', k"^) obtained from solving eq.(|T^. For comparison 
we also show gi{x, Q^) obtained from the standard LO QCD evolution (cf. eqs (|T2| , p!3[) ) as 
well as the input structure function gi'\x) given by eq. ([T9|) . We can see that the ln^(l/a;) 
resummation gives steeper structure function than that generated from the LO evolution 
and this effect is already visible for x < 10~^. 

The structure function (yfi(x, Q^) obtained from equations (^ and (|T7D still satisfies 
the Bjorken sum rule despite the fact that the double logarithmic terms might in principle 
lead to its violation. This can be seen as follows. One can show from eq.(^ that the first 
moment of the function gi is related to the first moment of /: 

i' "^'f^- = /„' ''^«!°'(^) + ir FfTTPA?) i' '''' 

It can easily be shown that the first moment of the function / vanishes and so the first 
moment of gi{x, Q"^) remains to be given by the first moment of gf'^ which is equal to gAl^-, 
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cf. eg. (PPD . Vanishing of the first moment of / follows from the fact that eq.(P!7|) would 
lead to the corresponding inhomogeneous integral equation for the moment function of / 
with the vanishing inhomogeneous term for the first moment. This inhomogeneous term 
is just the first moment of f^'^^ defined by eq. (|T8|) . The solution of such an equation then 
also vanishes. Vanishing of the first moment of f^^^ on the other hand follows from the 
fact that f^'^^ is proportional to the convolution of the splitting function Pgg with g^^\ cf. 
eq.(|l^), and the first moment of Pgg is equal to zero. 

Equations (H) and ([T7|) generate dynamically the small x extrapolation of gi{x,Q'^) 
and so one may in particular determine the contribution to the Bjorken integral coming 
from the region of very small values of x which is not being experimentally probed at 
present. We have thus estimated a contribution 

AI{xa, Xb, Q^) = r dxgiix, Q^) (22) 

to the Bjorken integral and found that A/(0, 0.003, 10 GeV^) = 0.0057 for the pure 
Altarelli-Parisi evolution calculations and 0.0080 for these with double logarithmic terms 
included. These numbers have to be compared with 0.004 obtained when a. gi = const 
behaviour, consistent with Regge prediction has been assumed and fitted to the lowest 
x data points for proton and deuteron targets (cf. An analogous number in the 

present case, A/(0, 0.003)=0.0026, where we have used gi*^ instead of gi in the integrand 
of A J, eq.(^. 
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Figure 2. Non-singlet part of the proton spin structure function gi(x, Q^). Points show the recent SMC 
results, [|[|3|, where g"'' = g\ — gi = 2{g\ — gf/{l — was obtained from the proton and deuteron 

data {ojd denotes a probability of the D state of the deuteron; ujd=0.05). Error bars are statistical; 
systematic errors are of the order of 20% of the statistical ones (in the lowest x bin) . Values of change 
with X as marked in the figure. Continuous line corresponds to the model calculations, eqs (|^), ( p^ ) and 
(19), performed at the measured values of {x,Q'^); broken line is a pure leading order Altarelli-Parisi 
prediction, eqs (|l^), ( [l^ ) and (^, at the same points. 

Let us finally comment upon the possible low extrapolation of the structure function 
gi{x,Q'^). For (for fixed pq) the structure function gi{x,Q'^) should be a finite 
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function of pg, free from any kinematical singularities or zeros at = 0. It may be 
seen that the structure function gi{x,Q'^) defined by eqs and with gi^\x) given 

by parametrisation (O) fulfills those criteria. This may not always hold for arbitrary 
parametrisation of gi which may include kinematical singularities at a; = 0. In 



such a case one may just use the same prescription as that which was adopted in ref. |24 
i.e. to replace in eq.(|) the function gf^\x) by gi^\x) where 



x = x[l + -^] (23) 




and to leave remaining parts of the calculation unchanged. After this simple rearrang- 
ment the structure function gi{x, Q^) can be extrapolated to the low region (for fixed 
2pq = Q'^/x) including the point = 0. It has however to be remembered that the 
(extrapolated) partonic contribution to the low region may not be the only one there. 




Figure 3. Non-singlet part of the proton spin structure function gi{x, Q"^) as a function of x for —10 
GeV^. Continuous line corresponds to the model calculations, eqs (^), ( p7| ) and (p^, broken line is a 
pure leading order Altarelli-Parisi prediction, eqs (^, (^^ and ([l9|), and a dotted one corresponds to 
9^1 \ eq-di). 



5. Summary and conclusions 

In this paper we have formulated the unified equation for the spin dependent non-singlet 
structure function g^'^{x,Q'^) which incorporated the LO Altarelli-Parisi evolution and 
the double logarithmic ln^(l/x) effects at small x. The equation treats on equal footing 
two potentially "large" logarithms ln{Q^) and ln(l/x). It also allows to extrapolate gi^ 
dynamically to the region of small x. 
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The non-singlet case has been chosen for its simphcity; dominant part of the double 
logarithmic ln^(l/x) terms and the QCD evolution here are generated by ladder diagrams 
with quark and antiquark exchange. It has also been chosen in order to meet the experi- 
mental need to know the low x behaviour of necessary when evaluating the Bjorken 
sum rule from the data. Elaboration of similar formalism for the singlet case is in progress. 

We solved that equation starting from the simple parametrisation of the non-perturbative 
part of gi^{x,Q'^). We found that the ln^(l/a;) effects are very significant already for 
X ^ 10"^ which may be probed in the possible future HERA measurements. We have also 
estimated the contribution from the small x region {x < 0.003) to the Bjorken sum rule 
and found it to be around 4% of the value of the sum. We have pointed out that the 
structure function g'^^{x,Q^) obtained from our formalism can easily be extrapolated to 
the region of low although it may not be the only contribution in this region. 
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